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We show that the boost-invariant and cylindrically asymmetric hydro- 
dynamic equations for baryon-free matter may be rewritten as only two 
coupled partial differential equations. In the case where the system ex- 
hibits the cross-over phase transition, the standard numerical methods may 
be applied to solve these equations. An example of our results describing 
non-central gold on gold collisions at RHIC energies is presented. 
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1. Introduction 

The success of the relativistic hydrodynamics in describing the RHIC 
data ^ HI El IU El El 13 suggests that the hot and dense matter produced at 
RHIC behaves like an almost perfect fluid [Sj- In this paper we discuss how 
the boost-invariant and cylindrically asymmetric relativistic hydrodynamic 
equations for baryon-free perfect fluid may be very conveniently reduced to 
only two coupled partial differential equations. Our presentation is based 
on the recent investigations published in Refs. [HI Ell- We argue that the 
effects of the cross-over phase transition may be included by the use of the 
temperature dependent sound velocity c s (T). As long as the function c s (T) 
satisfies the stability condition against the shock formation, the resulting 
equations may be solved with the help of standard numerical methods and 
used to describe the expansion of matter produced in the central region of 
ultra-relativistic heavy-ion collisions, such as studied in the present RHIC or 
future LHC experiments. The presented formalism is a direct generalization 
of the approach introduced by Baym, Friman, Blaizot, Soyeur, and Czyz 
jllj where the boost-invariant and cylindrically symmetric systems were 
considered, and the numerical calculations were performed only for the case 
of the constant sound velocity. 
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2. Hydrodynamic equations 

The relativistic hydrodynamic equations of the perfect fluid follow from 
the energy-momentum conservation law and the assumption of the local 
thermal equilibrium. For baryon-free matter they have the following form 

u^(Tu v ) = d»T, d^su") = 0, (1) 

where T is the temperature, s is the entropy density, and = 7(1, v) is 
the four-velocity of the fluid. For boost-invariant systems Eqs. Q may be 
reduced to three independent equations 

3d d 

— (rtsj) + — (rtsjv cos a) + — (ts'yv sin a) = 0, (2a) 
at Or ocp 

Q d d 

— (rT'yv) + r cos a— (T7) + sin a— (T7) = 0, (2b) 
at or ocp 

m 9 / da v sin a \ dT cos a dT , . 

T ^ v U + — " ma » + — w = a (2c) 



In Eqs. (|2*aj) - (|2*c|) the quantities t, r=y/x 2 + y 2 , and = tan l (y/x) are 
space-time coordinates which parameterize the plane z = 0. The quantity u 

is the magnitude of the fluid velocity, 7 = (l — v 2 ) 1 is the Lorentz factor, 
and a is the function describing direction of the flow, a = tan -1 (vt/vr) 

A 
dt 



see Fig. ^ The differential operator A in (|2c|) is defined by the expression 



d d d vsma d . . 

dt = di +VCOSa d-r + —W (3) 

Equations © are three differential equations for four unknown functions: 
T, s, v, and a. To solve them one has to use also the equation of state, i.e., 
the relation connecting T and s. In our approach, the equation of state is 
taken into account by the use of the temperature dependent sound velocity 

s{ ' 8e T 8s' [> 
and by the use of the potential $(T) defined by the differential equation 

d$ = — Cs dlns. (5) 

The form of the function c 2 (T) used in our calculations is defined and dis- 
cussed in more detail in Ref. [S]. Here we only note that the form of 
our sound velocity function automatically satisfies the condition against the 
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Fig. 1. Decomposition of the flow velocity vector in the plane 2 = 0. In our 
approach we use the magnitude of the flow v and the angle a as two independent 
quantities. 



shock formation [111 I14j . The integration of Eq. (JSJ allows us to express 
<3? in terms of the temperature, <3? = $y(T), or to express temperature in 
terms of T = T<j>(<]?). These two functions may be used to rewrite the 
hydrodynamic equations in the very concise form. In fact, with the help of 
the substitutions: 



tanh#, a± = exp(<l> ± 0), 



(6) 

the sum and the difference of the equations ()2a|) and ()2b|) may be rewritten 
as 

da± (v±c s ) da± (v±c s ) sinada± 

dt (1 ± c s v) dr (1 ± c s v) r dej) 

da cos a da N 



c s v 

(1 ± c s v) \ dr 

1 v cos a 

t r 



(1 ± c s v 

while the third equation in (|2cjl is 



a± = 0, 



a± 



(7) 



da 



(1 ~ v 2 ) 



sma 



d$ 

dr 



cos a 



<9$ 



9a sin a 9a sin a 
-u | cos a— H — H 



(8) 



Equations Q and © are three equations for three unknown functions: 
a_|_(i, r, <f>), a-(t,r,<fi), and a(t,r,<f>). We note, that the velocity i> and the 
potential $ are functions of a+ and a_, 

v = a+ ~ a " , $ = ^m(a + a_). (9) 
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a(r,4M)/ 


, a + (r,0,t) = a(r,^,t) 
\a_(r,(M) = a(-r,0,t) 
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Fig. 2. Construction of the functions a+(t, r, </>) and a_(i, r, 0) with the help of the 
function a(t,r,(f>). The function a(t,r,(j)) is symmetrically extended to negative 
values of r. 



Also the sound velocity may be expressed as the function of a + and a_, 
however, through a more complicated formula, see Ref. |10| . 

3. Boundary conditions 

In this paper we consider the collisions of identical nuclei with atomic 
number A which collide moving initially along the z-axis. The impact pa- 
rameter b points in the x-direction. The center of the first nucleus in the 
transverse plane is placed at xi = (xi,yi) = (—6/2,0), and of the second 
nucleus at X2 = (x2, 2/2) = (6/2, 0). In such a case we require that the mag- 
nitude of the flow vanishes at x = (x,y) = (0,0). This condition is fulfilled 
naturally by the ansatz 

a + (t,r,(f)) = a(t,r,(f>), r > 0, (10a) 
a-(t,r,4>) = a(t,-r,<j>), r < 0. (10b) 

We supplement the ansatz (|1U|) by the definition of the function a(t,r,cp) 
for the negative arguments, see Fig. [2J 

a{t,-r,4>) = a(t,r,4>), r > 0. (11) 

With the help of the definitions (|TU|) and (|lljl. Eqs. (JZJ) may be reduced to 
a single equation for the function a(t,r,(f>), 

da (v + c s ) da (v + c s ) sin a da 

at (1 + c s v) or (1 + c s v) r ocp 



c s v ( . da cos a da\ 

sin a— —— a + 



(1 + c s v) V dr r d(f> J (1 + c s v) [t 



1 v cos a 
- + 



a = 0, 



Here, similarly to the cylindrically symmetric case, the range of the variable 
r is extended to negative values. Eq. (jT3)) should be solved together with 
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Eq. l(HJl. where the range of r may be also extended trivially to negative 
values (definitions (|1U|) and imply that this equation is in fact invariant 
under transformation: r — > —r). The use of the polar coordinates in the 
transverse plane requires also that all functions at (f> = and <j> = 2ir are 
equal: a(i, r, 0) = a(t, r, 2ir), a(t, r, 0) = a(i, r, 2ir). 



4. Initial conditions 

In the following we assume that the hydrodynamic evolution starts at 
a typical time t = to = 1 fm. We assume also that the initial temperature 
profile is connected with the number of participating nucleons 

/ \ V 3 

T(t ,x,y) = const . (13) 

The idea to use Eq. (|13|) follows from the assumption that the initially 
produced entropy density a(to,x,y) is proportional to the number of the 
nucleons participating in the collision at the position (x,y). Since the con- 
sidered systems are initially very hot (with the temperature exceeding the 
critical temperature T c ), they may be considered as systems of massless 
particles, where the entropy density is proportional to the third power of 
the temperature. In this way we arrive at Eq. (|13j) . The initial form of the 
functions v(t,r,cj)) and a(t,r, <f>) is 

H T 

v(t ,r,(j)) = v (r) = a(t ,r,4>) =0, (14) 

<1 + Hlr 2 



where Hq is a parameter defining the initial transverse flow formed in the 
pre-equilibrium stage. In the present calculations we use a very small value 
Hq = 0.001 fm . The two initial conditions, Eqs. (|13fl and (JHJ), m ay be 
included in the initial form of the function a(r) if we define 



a(t = t ,r,<f>) = exp < $ T 



const 



( dN p \ 
\dxdy J 



Vl-u"(r)' 



5. Results 

In this section we give an example of our results describing the collisions 
with the impact parameter b = 7.6 fm, typical for centrality class c = 0—80% 
[El 1 . In this case the initial central temperature To equals 2.5 T c , where 
T c is the critical temperature, Tq = T(to,0,0). The part a) of Fig. |31 shows 



1 In the paper [TH] we considered the centrality classes: 0-20%, 20-40%, 40-60% 
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the temperature profiles for different values of time: t = 1, 4, 7, 10, 13, 16 fm. 
The solid lines represent the temperature profiles along the x-axis (</> = 0), 
while the dashed lines represent the profiles along the y-axis (<j) = 7r/2). 
One can observe that during the whole considered evolution time, 1 fm 
< t < 16 fm, the y-extension of the system remains larger than the x- 
extension. However, the relative magnitude of this effect decreases with 
time, indicating that the cylindrical symmetry is gradually restored as the 
time increases. In this part of Fig. |3] one can also notice the effect of the 
phase transition; initially the system cools down rather rapidly, later the 
cooling down is delayed and the main visible effect is the increase of the 
volume of the system. We note that this behavior is related to the sudden 
decrease of the sound velocity in the region T ~ T c . The part b) of Fig. 
121 shows the isotherms in the t — r space, again for (f> = (solid lines) and 
(j) = it/2 (dashed lines). The pairs of isotherms indicate different values of 
the temperature. They start at T = 1.8 T c and go down to T = 0.2 T c , with 
a step of 0.2 T c . It is interesting to observe that the solid and dashed lines 
cross each other. This effect means that the central (relatively hotter) part 
of the system acquires a pumpkin-like shape during the evolution of the 
system. Such pumpkin-like regions, however, shrink and disappear during 
further expansion. In the part c) of Fig. |3] the profiles of the function 
a(t, r, (j)) are shown for t = 1,4, 7, 10, 13, 16 fm and <p = 7r/4. For = and 
<fi = vr/2 the function a(t,r, 4>) vanishes due to the symmetry reasons. In the 
first and third quadrant the values of a are predominantly negative, while in 
the second and fourth quadrant the values of a are positive. This behavior 
characterizes the direction of the flow which has the tendency to change the 
initial almond shape into a cylindrically symmetric shape. In the part d) 
the velocity profiles are shown, again for t = 1, 4, 10, 16 fm. Similarly to 
the previously discussed figures, the solid lines are the profiles for (j) = 
(x-direction), whereas the dashed lines are the profiles for (p = ir/2 {y- 
direction). One can observe that the magnitude in the x-direction is larger 
than the magnitude in the y-direction, which is an expected hydrodynamic 
behavior caused by larger pressure gradients in the x-direction. Exactly this 
effect is responsible for the observed azimuthal asymmetry of the transverse- 
momentum spectra, quantified by the values of the t>2 coefficient. Finally, in 
the part e) we show the contour lines of the temperature, again for different 
values of time. These plots visualize the time development of the system. 
The arrows describe the magnitude and the direction of the flow (for better 
recognition the angle a is magnified by a factor of 3). 
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6. Summary 

In this paper we presented a new and concise treatment of the boost- 
invariant and cylindrically asymmetric relativistic hydrodynamic equations. 
The presented formalism is a direct generalization of the approach intro- 
duced by Baym et al. in Ref. In the studied case, the symmetry 
of the problem allows us to rewrite the hydrodynamic equations as only 
two coupled partial differential equations, (jHJ) and (jlMjl . which automati- 
cally lead to the fulfillment of the requested boundary conditions for the 
velocity and the temperature at the center of the system. The effects of 
the phase transition are included in this scheme by the use of the tempera- 
ture dependent sound velocity. The presented results of the hydrodynamic 
calculations, supplemented with the appropriate freeze-out model (e.g., the 
single-freeze-out model of Refs. HS1E1), ma y ^ e used to calculate different 
physical observables. To achieve this task, in the closest future we intend 
to combine our hydrodynamic approach with the statistical Monte-Carlo 
model THERMINATOR [Tgj . 
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Fig. 3. Time development of matter characterized by the initial conditions (|14|) - 
(Q2J with H a = 0.001im~\ T = T(t ,0,0) = 2.5 T c , and b = 7.6 fm. 



